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We model analytically the dynamics of a cytoskeletal filament in a motility assay. The filament 
is described as rigid rod free to slide in two dimensions. The motor proteins consist of polymeric 
tails tethered to the plane and modeled as linear springs and motor heads that bind to the filament. 
As in related models of rigid and soft two-state motors, the binding/unbinding dynamics of the 
motor heads and the dependence of the transition rates on the load exerted by the motor tails play 
a crucial role in controlling the filament's dynamics. Our work shows that the filament effectively 
behaves as a self-propelled rod at long times, but with non-Markovian noise sources arising from the 
coupling to the motor binding/unbinding dynamics. The effective propulsion force of the filament 
and the active renormalization of the various friction and diffusion constants are calculated in terms 
of microscopic motor and filament parameters. These quantities could be probed by optical force 
microscopy. 



There has recently been renewed interest in motility as- 
says where semifLexible actin filaments are driven to slide 
over a "bed" of myosin molecular motors. Recent exper- 
iments at high actin density have revealed that the col- 
lective behavior of this simple active system is very rich, 
with propagating density waves and large scale-swirling 
motion [HE], not unlike those observed in dense bacterial 
suspensions [3]. In an actin motility assay the polymeric 
tails of myosin motor proteins are anchored to a surface, 
while their heads can bind to actin filaments jl]. Once 
bound, the motor head exerts forces and drives the fila- 
ment's motion. This system provides possibly the sim- 
plest realization of an active system that allows detailed 
semi-microscopic modeling. 

Stochastic models of the collective action of motor pro- 
teins on cytoskeletal filaments in one dimension have 
been considered before by several authors, with empha- 
sis on the acto-myosin system in muscles and on the mi- 
totic spindle [S]. When working against an elastic load, 
the motor assemblies have been shown to drive periodic 
spontaneous activity in the form of oscillatory instabil- 
ities, which in turn have been observed ubiquitously in 
a variety of biological systems [BHTU] . These instabilities 
arise in the model from the collective action of the mo- 
tors and the breaking of detailed balance in their dynam- 
ics and manifest themselves as a negative effective fric- 
tion of the filament. When free to slide under the action 
of an external force, the filament can exhibit bistability 
that manifests itself as hysteresis in the force velocity- 
curve [TTJ HI] . A large body of earlier work has mod- 
eled the motors as rigid two-state systems attached to 
a backbone and bound by the periodic potential exerted 
by the filament on the motor head [6l [11] [13] . In a sec- 
ond class of models the motors have been modeled as 
fiexible springs [TU |TS| . The motor heads bind to the fil- 
ament and unbind at a load-dependent rate. In this case 
the dynamic instability arises from the dependence of the 
unbinding rate on the tension exerted by springs |16H18j . 



Recent work by Guerin et al. TS] has generalized the two- 
state model by taking into account the fiexibility of the 
motors, showing that both models can be obtained in a 
unified manner for different values of a parameter that 
compares the stiffness of the motors to the stiffness of the 
periodic potential provided by the filament. 

In this paper we consider a model of a rigid filament 
free to slide in two dimensions under the driving action of 
motor proteins uniformly tethered to a two-dimensional 
plane. The model considered is a modification of the 
"crossbridge" model first introduced by Huxley in 1957 
to describe motor-induced contractile behavior of muscle 
fibers [3D|. The motor proteins' polymeric tails are mod- 
eled as linear springs that pull back on the bound motor 
heads. After attachment, the motor heads slide along the 
filament at a velocity that depends on the load exerted by 
the fiexible motor tails. The sliding and subsequent de- 
tachment play the role of the motor's power stroke. The 
binding/unbinding dynamics of the motor heads and the 
dependence of the transition rates on the load exerted 
by the motor tails play a crucial role in controlling the 
dynamics of the fie, effectively yielding non-Markovian 
noise sources on the filament. Related models have been 
studied numerically [HI [THI HI] ■ The results presented 
here are obtained by generalizing to two dimensions the 
mean field approximation for the motor dynamics de- 
scribed for instance in Ref. The mean-field theory 
neglects convective nonlinearities in the equation for the 
probability of bound motors and correlations in the mo- 
tors on/off dynamics, but it is expected to be adequate 
on time scales large compared to that of the motor on/off 
dynamics and for a large number of motors. This is sup- 
ported by the results of [11^ for a model of rigid two-state 
motors. 

We begin by revisiting the one-dimensional problem. 
We discuss the steady-state response of the filament to 
an external force and present new results on the dynam- 
ics of fiuctuations about the sliding steady state. The 
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force- velocity curve is evaluated analytically and exhibits 
bistability and hysteresis, as obtained in Ref. ^HJ for a 
rigid two-state motor model. A new result is an expres- 
sion for the effective propulsion force on the filament due 
to the motors in terms of physical parameters character- 
izing the motor proteins. Next, we analyze the fluctu- 
ations about the steady state by evaluating the mean- 
square displacement of the filament. We show that the 
coupling to the motor binding/unbinding dynamics yields 
non-Markovian noise sources with time correlations con- 
trolled by the duration of the motors' binding/unbindig 
cycle. Since the filament has a finite motor-induced ve- 
locity even in the absence of applied force, the mean- 
square displacement is ballistic at long time. The fiuctu- 
ations of displacement about this sliding state are, how- 
ever, diffusive at long times with an enhanced diffusion 
constant. This enhancement is controlled by the depen- 
dence of the motors' unbinding rate on the load exerted 
on the bound motors' heads by the tethered tails and 
vanishes for unloaded motors. 

We then consider the case of a filament in two dimen- 
sions, to analyze the effect of the coupling of transla- 
tional and rotational degrees of freedom in controlling 
the dynamics. At steady state, motors yield an effective 
propulsion force along the long axis of the filament, as in 
one dimension, but no effective torque. This is in con- 
trast to phenomenological models considered in the lit- 
erature [55] that have considered the dynamics of active 
rod-like particles in the presence of both effective inter- 
nal forces and torques. As a result, in the steady-state 
the filament slides along its long axis and the dynam- 
ics in this direction is essentially one dimensional, with 
a motor-induced negative friction instability and bista- 
bility and hysteresis in the response to an external force. 
Motors do enhance both the transverse and the rotational 
friction coefficients of the filament. The enhancement 
of rotational friction could be probed by measuring the 
response to an external torque. Since the finite motor- 
induced propulsion is along the filament axis, whose di- 
rection is in turn randomized by rotational diffusion, the 
mean velocity of the filament is zero in the absence of 
external force, unlike in the one-dimensional case. The 
mean square displacement is therefore diffusive at long 
times, with behavior controlled by the interplay of non- 
Markovian effects due to the coupling to motor dynamics 
with coupled translational and rotational diffusions. The 
filament performs a persistent random walk that consists 
of ballistic excursions at the motor-induced propulsion 
speed, randomized by both rotational diffusion and the 
motor binding/undinding dynamics. The crossover to 
the long-time diffusive behavior is controlled by the inter- 
play of motor-renormalized diffusion rate and duration of 
the motor binding/unbinding cycle. The effective diffu- 
sion constant is calculated in terms of microscopic motor 
and filament parameters. Its dependence on activity, as 
characterized by the rate of ATP consumption, could be 
probed in actin assays. 

Finally, our work provides a microscopic justification 



of a simple model used in the literature that de- 
scribes a cytoskeletal filament interacting with motor 
proteins tethered to a plane as a "self-propelled" rod, 
although it also shows that the effective noise is ren- 
dered non-Markovian by the coupling to the motors' 
binding/unbing dynamics. It also provides microscopic 
expressions for the self-propulsion force and the various 
friction coefficients in terms of motor and filament pa- 
rameters and shows that this effective model fails beyond 
a critical value of motor activity, where the effective fric- 
tion changes sign and the filament exhibits bistability and 
hysteresis. 

I. THE MODEL 

In our model the motor proteins are described as com- 
posed of polymeric tails attached permanently to a two- 
dimensional fixed substrate and motor heads that can 
bind reversibly to the filament. Once bound, a motor 
head moves along the filament thereby stretching the tail. 
This gives rise to a load force on the motor head and on 
the filament. Eventually excessive load leads to detach- 
ment of the motor head. 



A. Filament dynamics 

The actin filament is modeled as a rigid polar rod of 
length L that can slide in two dimensions. It is described 
by the position r of its center of mass and a unit vector 
u = (cos(6'), sin(0)) directed along the rod's long axis 
away from the polar direction of the rod, which is in 
turn defined as the direction of motion of bound motors. 
In other words, bound motors move along the rod in the 
direction — u. In contrast to most previous work [6l [3 
1131 119j , and given our interest in modeling actin motility 
assays, we assume the substrate is fixed and consider the 
dynamics of the filament. Our goal is to understand the 
role of the cooperative driving by motors in controlling 
the coupled rotational and translational dynamics of the 
rod. 

The dynamics of the filament is described by coupled 
equations for the translational and orientational degrees 
of freedom, given by 

£ • 9tr = Fa + Fe,t + il{t) , (la) 
^edtO = T, + T,^t+Tle[t) . (lb) 

Here we have grouped the forces and torques into the 
effects due to the motors, i.e. the activity. Fa and Ta, 
external forces and torques Fext an Text and the stochas- 
tic noise not due to motors. The friction tensor is given 
by ^ = C||uu -I- Cj_ (5 — uu) with Cy and C_l the fric- 
tion coefficients for motion longitudinal and transverse 
to the long direction of the rod, and Ce is the rotational 
friction coefficient. For the case of a long, thin rod of 
interest here, = C,^_/2. The random force ri{t) and 
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FIG. 1: The figure shows the four steps of a motor cycle. In 

(a) a filament is sliding with velocity v over a uniform density 
of unbound motors with tails tethered to the substrate. In 

(b) a motor attaches to the filament at a position so from the 
filament's mid-point. The stretch of the motor tails at the 
time of attachment is neglected. In (c) the motor has walked 
towards the polar head of the filament, stretching the tails by 
an amount A. Finally, in (d) the bound motor detaches and 
relaxes instantaneously to its unstretched state. The filament 
has undergone a net displacement in the direction opposite to 
that of motor motion. 

random torque T]e{t) describe noise in the system, in- 
cluding nonthermal noise sources. For simplicity we as- 
sume that both r}{t) and rjoit) describe Gaussian white 
noise, with zero mean and correlations {rji{t)rjj{t')) — 
2B,jS{t - t') and {ri0{t)r]g{t')) = 2BeS{t - t% where 
Bij = B\\UiUj + Bj^ ((5,y - iiiUi). 



B. Individual motor dynamics 

We model the interaction cycle of an individual motor 
protein with the filament as shown in Fig. [l] for a one- 
dimensional system. The tail of a specific motor is fixed 
at position xt in the plane. At a time to the head of this 
motor attaches to a point on the filament. The position 
of the motor head at the time of attachment is x/i(io) — 
r(io) +sou(to): where r(io) and u{to) denote the position 
of the center of the filament and its orientation t = to 



and So G ["-'^/^j^/S] parametrizes the distance of the 
point of attachment from the center of the filament (cf. 
Fig. [TJb)). We assume that motor proteins will attach to 
parts of the filament which are within a distance of the 
order of the size of the motor protein. The stretch of the 
motor tail at the time of attachment is then of order of 
the motor size and will be neglected, i.e. x^(to) — x^ = 
0, or motors attach to the part of the filament directly 
overhead without any initial stretch. 

For t > to the motor head remains attached to the fil- 
ament and walks along it towards the polar head (— u di- 
rection) until detachment. The tails, modeled as a linear 
spring of force constant k, exert a load f — —kA{t, r; sq) 
on the head, where A(i, t; sq) — Xh{t) — Xt is the stretch 
at time t of a, motor protein that has been attached for 
a time r, i.e. t = to + t (cf. Fig. [ijc)). Since we assume 
A(to) = 0, we can also write 

A(i,T;so) = r{t)-r{t-T) + a{t,T)u{t) 

+.S0 [u(<) - uit - r)] , (2) 

where a{t, r) — s{t) — s{t — t) is the distance traveled 
along the filament at time t by a motor head that has 
been attached for a time r, measured from the initial at- 
tachment position. So- The kinematic constraint imposed 
by the condition of attachment requires 

dt^{t, r; So) - v(<) - v(t - r) + u(<) [vra{t) - - r)] 

+n{t)a{t,T) + so[n{t)-n{t-T)] , (a) 

where n{t) = dtu.{t) = 6n{t) is the angular velocity of 
the rod and Vm{t) = dts{t) the velocity of the motor head 
along the filament. We have introduced a unit vector 
fi = z X u normal to the long axis of the filament. Then 
(z,u, n) defines a right-handed coordinate system with 
in-plane axes longitudinal and transverse to the filament. 
We note that Eq. ([S]) can also be written as 

dtA{t,T;so) + drA{t,T;so)^v{t)+v^{t)u{t) 

+n{t)a{t,T) + .son{t) . (4) 

While the motor remains bound, the dynamics of the 
motor head along the filament is described by an over- 
damped equation of motion 

um = + u • f (5) 

where > is the stall force, defined as the force where 
the velocity Vm = s of the loaded motor vanishes. Since 
motors move in the — u direction, generally Vm = s < 0. 
Letting /|| = u • f = — fcA|[, Eq. ([s]) can also be written 
as 

Vm[t) = -Vo (1 j ' 

where vo = /s/Cm ^ A/i > is the load-free stepping ve- 
locity, with A/1 the rate of ATP consumption. The motor 
velocity is shown in Fig. (pi) as a function of the load /n . 
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FIG. 2: The velocity —Vm of a loaded motor head as a function 
of the load /y = u • A. The figure shows the stall force fa 
where Vm ~ and the detachment force —fd- 

The motor head velocity also vanishes for /y < —fd, when 
the motor detaches. The linear force- velocity relation for 
an individual motor is consistent vi^ith experiments on 
single kinesin molecules [24] . 

The active force and torque on the filament due to an 
individual bound motor can then be expressed in terms 
of these quantities as 

ia{t,T;sa) ^ -kA{t,T;so) , (7a) 

Ta{t, t; So) = -z • [(so + cr(i, T))u{t) X kA{t, r; so)] . 

(7b) 

Finally, after traveling along the filament for a time 
"'detach I the motor head detaches and the head position 
relaxes instantaneously back to the fixed position xt of 
the tail. 

We note that we shall not be considering the possibil- 
ity of direct interactions of motors with each other. We 
have also not considered stochastic aspects of the motor 
motion along the filament (Eq. ([s])). 

C. Motor binding and unbinding 

Next we need to describe the stochastic bind- 
ing/unbinding dynamics of the motor heads. We assume 
the motor tails are attached to the substrate with a ho- 
mogeneous surface density pm, such that for a rod of 
length L and width b a maximum of iV = PmLb motors 
can be bound at any given time. Following Guerin et 
al. [IH], we denote by Vb{t,T;so) the probability that a 
motor head that has attached at so at a time to, has 
remained attached for a duration r at time t. For sim- 
plicity in the following we assume that the probability 
that a motor attaches at any point along the filament is 
uniform, i.e., Vb{t, r; sq) ~ j^Pb{t, r). We further assume 
that when motors unbind they relax instantaneously to 
the unstretched state. The time evolution of the binding 



probability is then given by 

dtPbit,T) + drPb{t,T) = - {uJu{A{T))),„Pb{t,T) 

+ UJb5iT)puit) , (8) 

where Pu{t) is the probability that a motor be unbound 
at time t. The probability distribution is normalized ac- 
cording to 

/•oo rL/2 

/ dT dsoVb(t,T;so)+pu{t)^l . (9) 

JQ J-L/2 

In Eq. ([s]), uJu{A{t)) and ojb are the rates at which a mo- 
tor head with tails stretched by an amount A(t,T) un- 
binds from and binds to the filament, respectively. The 
binding rate LUb will be assumed to be constant. In con- 
trast, the unbinding rate uiu is a strong function of the 
stretch of the motor tails, that has to be obtained by 
solving Eq. Q, with initial condition A(t = 0,t) = 0. 
We will see below that the nonlinear dependence of 
the unbinding rate plays an important role in control- 
ling the filament dynamics. In two dimensions the un- 
binding rate lOu also depends on the initial attachment 
point So along the filament. To be consistent with our 
ansatz that the probability that the motor attaches at 
any point along the filament is uniform, we have replaced 
the rate in Eq. Q with its mean value (w„)so, where 
(...)so = /-L/2 if"- denotes an average over the initial 
attachment points. 

The unbinding rate is controlled by the work done by 
the force (load) acting on the motor head, which in turn 
is a linear function of the stretch A. A form that has 
been used extensively in the literature for one-dimensinal 
models is an exponential, a;„ = woe"'^', where u>o is the 
unbinding rate of an unloaded motor and a is a charac- 
teristic length scale that control the maximum stretch of 
the tails above which the motor unbinds 34 . The expo- 
nential form represents an approximation for the result 
of a detailed calculation of the average time that a mo- 
tor moving along a polar filament spends attached to the 
filament as a function of a tangentially applied load [25] 
and is consistent with experiments on kinesin [26] . This 
form can easily be generalized to to the case of a fil- 
ament sliding in two dimensions where the motor load 
had both components tangential and transverse to the 
filament. It is, however, shown in the Appendix that 
within the mean-field approximation used below the ex- 
ponential form yields a steady-state stretch A that satu- 
rates to a finite value at large velocity v of the filament. 
This is unphysical as it does not incorporate the cutoff 
described by the detachment force in Fig. [2] For this 
reason in the mean-field treatment described below we 
use a parabolic form for the unbinding rate as a function 
of stretch, 

c^„(A)=a.o[l + «'|Ap] , (10) 

where for simplicity we have assumed an isotropic de- 
pendence on the magnitude of the stretch in terms of a 
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single length scale, a~^. An explicit comparison of the 
two expressions for the unbinding rates is given in the 
Appendix. 



II. MEAN FIELD APPROXIMATION 

To proceed, we introduce several approximations for 
the motor dynamics. First, we restrict ourselves to the 
dynamics on times scales large compared to the attach- 
ment time T of individual motors. For t 3> r we approx- 
imate 

a{t,T)^v„,{t)T , (12a) 
Ait, r; So) ^ [v(t) + i;™(<)u(i) + son{t)] r . (12b) 

This approximation becomes exact for steady states 
where the filament and motor velocities are independent 
of time. We also stress that in Eqs. ( |12a[ ) and ( |12b[ ) a and 
A are still nonlinear functions of t due to the dependence 
of Vm on the load force. 

Secondly, we recall that we have assumed that the at- 
tachment positions Sq are uniformly distributed along the 
filament and can be treated as independent of the resi- 
dence times r. Finally, we make a mean field assump- 
tion on the probability distribution of attachment times, 
which is chosen of the form P{t,T) = S{t — TMF)Pb{t), 
with Pb{t) the probability that a motor be attached at 



In the following we will work in the mean-field approxi- 
mation and remove the label MF from the various quan- 
tities. 



III. ACTIVE FILAMENT SLIDING IN ONE 
DIMENSION 

We first consider the simplest theoretical realization of 
a motility assay experiment, where the actin filament is 
sliding over a one dimensional track of tethered motor 



The total active force and torque on the filament aver- 
aged over the original positions and the times of attach- 
ment can be written as 

(11a) 
(lib) 
I 

time t regardless of the its attachment time. The mean- 
field value of the attachment time is determined by re- 
quiring 

tmf = (A(rMF)))so]~^ • (13) 

In previous literature a similar mean field assumption 
has been stated in terms of the stretch, A |71 [S]. In 
the present problem, however, where filaments can slide 
in two dimensions, it is necessary to restate the mean- 
field theory in terms of the residence time t as the active 
forces and torques depend on both the stretch A of the 
motor tails and the distance a traveled by a bound motor 
head along the filament. These two quantities are in turn 
both controlled by a single stochastic variable, identified 
with the residence time r. The rate of change of the 
probability Pb(t) that a motor be bound at time t is then 
described by the equation 

dtPbit) - -r^lpbit) Pb{t)] , (14) 

The mean field active force and torque due to the mo- 
tors are then given by 



(15) 
(16) 

I 

proteins. A closely related model, where the filament is 
elastically coupled to a network, has been used exten- 
sively in the literature to study the onset of spontaneous 
oscillations arising from the collective action of the bound 
motors [51 [71 [TT] . Previous studies of freely sliding fila- 
ments, as appropriate for the modeling of motility assays, 
have also been carried out both analytically and numer- 
ically |18j . Our work contains some new results on the 
response to an external force of a filament free to slide 
under the action of active crosslinkers and also on the 
filament fluctuations. 



/>oo 

Fait) = -Nk / dT {Pbit, t) A{t, r; so)).„ , 
Jq 

Ta{t) = -Nk / dT {Pb{t, r) z • [(so + a{t, r))u(t) x A(t, r; so)]),„ 



Ff'' ^ -kN{A{t,TMF;sMt)},„ , 

F = -kN{pb{t) z ■ [(so + a{t, TMF))u{t) X A(i, tmf; sq)]) 
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Parameters 


Myosin-II 


Kinesin 




~ 2 nm 


~ 8 nm 


5s 


~ 1 nm 


~ 25 nm 


e 


~ 2 


~ 0.32 



TABLE I: Typical values of the length scales £o ~ «o/i^o and 
5a ~ fs/k introduced in the text and the ratio e for myosin 
II and kinesin. The parameters are taken from Refs. [27] 
and [B]. 



The Langevin equation for the center of mass coordi- 
nate X of the filament is given by 



Ci = Fa(t)+Fe,t+r;(t) , 



(17) 



where x is the center-of-mass velocity of the filament and 
the mean-field active force is given by 

^^F(t) = -fciVpf,(<)A(i,T). 



In o ne di mension the dependence on sq drops out and 
Eq. (12b) simply gives A ~ (a; -|- Vm)T- Substituting 
Eq. (|6| for Vm, we can solve for A as a function of a; and 



A(i,T) = 



(i - vo)/uja 



(19) 



and Ec^. ( 13 ) for the mean attachment time becomes 



[f-i(i) + e]^ 



(20) 



where f — ujqt and e = kv^/ fgUjQ. The parameter e is the 
ratio of the length £o = vq/loq traveled by an unloaded 
motor that remains attached for a time to the stretch 

= fs/k of the motor tails at the stall force, fs- Typical 
values for these length scales and the parameter e are 
given in Table [ij 

It is convenient to rewrite the mean residence time f 

as 



1 



(21) 



where u — x/vq and we have introduced a dimensionless 
parameter v — ia that controls the dependence of the 
unbinding rate on the load exerted on the bound heads 
by the stretched motor tails, with 



1 

To 



1 



(22) 



the geometric mean of the two length scales introduced 
earlier. For stiff motors, with 1 or Ss, i^Sg, 

while for floppy, easy to stretch motors, corresponding to 
e <C 1 OY £q <^ Sg, £ ^ £(). Setting v = corresponds 
to neglecting the load dependence of the unbinding rate. 
The exact solution to Eq. (21) for the mean residence 



time f{x) as a function of the filament velocity can be 
determined and is discussed in the Appendix. Clearly r 
has a maximum value at x — vq, where r = u>q^ and 
decays rapidly as |i — uo| grows. 



A. Steady State and its Stability 

We begin by characterizing the steady state dynamics 
of the filament in the absence of noise. Incorporating for 
generality an external force F^xt, the steady state velocity 
V of the filament is obtained from the solution of the 
nonlinear equation 



(v ^ F,xt + Fa{v) 



(23) 



where i^a(u) = —kNpi,s{v)A{v). The steady state stretch 
A(w) is given by Eq. (19) with x — v and 



UJbT{v) 
1 + T{v)uih 



(24) 



with t{v) given by Eq. ( |21[ ) for x = v. To gain some 
insight in the behavior of the system, we expand the 
active force as Fa{v) ~ Fp + {^)^^qV + 0{v^), with 
Fp = Fa{v = 0). Retaining only terms linear in v this 
gives a steady state force/velocity relation of the form 



(C + Ca)v = Fe,t + Fp 

with a filament "propulsion" force Fp 

Npbsok£o 



F„ 



(25) 



(26) 



where pbso = ''/[^ + (l^^)^o with r = u}b/{iJ->o+^b) the 
duty ratio, and tq = f{v = 0). The active contribution 
Cq = - (^)„^o *° ^^'^ friction is given by 



fc|Ao| 



PbsO- 



l-r\ 2a'^Al£o 



where Aq = A(u = 0) — —<-D/yiQ 
of external force, the filament will slide at a velocity 



£0/(7^0 ^ + e)- In the absence 



Vs=Fp/{C + Ca) 



(28) 



due to the action of the motor proteins. This motion 
is in the polar direction of the filament and opposite to 
the direction of motion of bound motors along the fila- 
ment. Phenomenological models of motility assays have 
described the actin filaments as "self-propelled" Brown- 
ian rods. Our model yields a microscopic calculation of 
such a "self-propulsion" force Fp in terms of microscopic 
parameters characterizing the motor proteins. We note 
that —Fp can also be interpreted as the "stall force" of 
the filament, i.e. the value of Fext required to yield ?; = 0. 
This is a quantity that may be experimentally accessible 
using optical force microscopy. 

If we neglect the load dependence of the unbinding rate 
by letting j/ = 0, the mean number of bound motors is 
simply Nr and Fp = Nrk£, with £ given by Eq. ( 22 ) . In 
this limit the sliding velocity in the absence of external 
force can be written as 



vo 



1 + C/CS 



(29) 
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where the active friction C° = Nrki/vo > is always 
positive. The sUding velocity vanishes when — and 
it saturates to its largest value Vq when the number Nr 
of bound motors becomes very large and ^ C- The 
behavior is controlled by the parameter e. If the motors 
are easy to stretch, i.e., e ^ 1, then the propulsion force is 
determined entirely by the elastic forces exerted by these 
weak bound motors, with ~ Nrkio. On the other 
hand stiff motors, with e 3> 1, stall before detaching. 
The propulsion force is then controlled by the motor stall 
force, with F° ~ Nrf,. 

The load-dependence of the unbinding rate changes 
qualitatively the behavior of the system. In particular, 
the net friction C + Ca can become negative, rendering 
the steady state unstable. This instability was already 
noted in Ref. [11] for a two-state model of active linkers 
and in Ref. [12] for a two state "soft" motor model. The 
full nonlinear force- velocity curves are shown in Fig.[3]for 
various values of the motor stiffness k, for parameters ap- 
propriate for acto-myosin systems. In the steady state, as 
we increase the active parameter k while keeping the sub- 
strate friction constant, the Fcxt — v curve becomes non- 
monotonic, and two distinct regions of bistability emerge. 
To understand the increase of the bistability region with 
motor stiffness, we note that the active force is simply 
proportional to k, hence naively one would indeed expect 
its effect to be more pronounced for stiff motors. The de- 
tailed behavior is, however, determined by the interplay 
of the mean residence time r that motors spend bound 
to the filament and the stretch, A. Soft, floppy motors 
have large stretches, controlled mainly be the length £o 
traveled by an unloaded motors advancing at speed vq. 
On the other hand, their residence time is small and the 
overall effect of the active force remains small. In con- 
trast, stiff motors have a small stretch, of order of the 
stretch Sg = fs/k oi a stalled motor, but long residence 
times and are collectively capable of slowing down the 
filament and even holding it in place against the action 
of the external force, driving the negative friction insta- 
bility. At even larger values of the external force motors 
are effectively always unbound due to the fast sliding of 
the filament and the velocity-force curve approaches the 
linear form obtained when no motors are present. This 
behavior is best seen from Fig. [5] 

The region of non-monotonicity of the force-velocity 
curve and associated bistability can also be displayed as 
a phase diagram, as shown in Fig. [4] The stiffness of 
myosins is about 5 pN/nm and the actin filament friction 
was estimated to be of order 0.003 pNs/nm in Ref 4J. 
In actomyosin systems the negative friction instability 
should therefore be observable in a range of experimen- 
tally relevant parameters. Kinesin motors have floppier 
tails and a smaller stiffness of about 0.5 pN/nm. In this 
case bistability effects should be prevalent only at very 
low filament friction, C ^ 0.001 pNs/nm. A proper es- 
timate of the region of parameters where the instability 
may be observable is rendered difficult by the fact that 
the onset of negative friction is also a strong function of 
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10000 
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FIG. 3: (Color online) Force-velocity curves for ( = 
0.002 pNnm"^s and various values of the motor stiffness k, 
showing the transition to non-monotonicity as k increases. 
The values of the stiffness k (in pN/nm) and the correspond- 
ing values for (in nm) and e are as follows: A; = 0, 
a"^ = 0, e = (black dotted line); fc = 1 , a"^ = 0.75, e = 0.5 
(red dashed line); k = 2, = 1.5, e = 1 (blue dashed- 
dotted line); k — 8, = 6, e = 4 (black solid line). At high 
velocities the curves merge into the linear curve _Foxt = 
(black dotted line), corresponding to the case where no mo- 
tors are present. The remaining parameters have the follow- 
ing values: = pmLb = 100, vo = 1000 nm/s, /s = 4 pN, 



0.5 (ms) 



0.06. 



the density of motors tethered to the substrate, which in 
turn affects the value of the friction In general, we 
expect that a high motor density will be needed for the 
instability to occur. On the other hand, if the density of 
motors is too high, the friction C will be enhanced and 
the instability suppressed. 

We stress that the force-velocity curves displayed in 
Fig. [3] have been obtained by calculating Fext as a func- 
tion of V. In an experiment one would tune the applied 
force and measure the resulting velocity. The system 
would not access the unstable regions of negative friction, 
but rather follow the hysteretic path sketched in Fig. [5] 
The discontinuous jump may occur at the boundary of 
the stability region, as shown in the figure, or before such 
a boundary is reached, corresponding to what is known 
as "early switching" . 

To summarize, motors have two important effects on 
the steady state dynamics of the filament. First, they 
make the filament self-propelled, in the sense that in the 
absence of an external force the filament will slide at a 
velocity Vs given by Eq. ( 28 1 . The value of Vs increases 



with increasing motor stiffness and of course vanishes for 
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FIG. 4: (Color online) "Phase diagram" in k-C, plane show- 
ing the region where the F^^t-v curves exhibit non-monotonic 
behavior (blue shaded region) for N — pmLb — 100 and 
Wo = 1 s"'^, fs~4: pN, a/k — 1.33 pN, cjq — 0.5 (ms)""*^, 
r = 0.06. 

V (nm/s) 
4000 ^ 
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FIG. 5: The figure sketches the hysteretic behavior that may 
be obtained in an experiment where an external force F^^t 
is applied to a filament in a motility assay. The response of 
the filament will generally display two regions of hysteresis, 
at positive and negative forces. 



vq = 0, corresponding to the vanishing of the rate of 
ATP consumption Afi. The sliding velocity Vs is shown 
in Fig. [6] as a function of the parameter e inversely pro- 
portional to the motor stall force for a few values of the 
maximum number of motors that can bind to the fila- 
ment. A second important effect of motor activity is the 
discontinuous and hysteretic response to an external force 
displayed in Fig. [5j When Foxt = the filament slides at 
the motor- induced velocity Vg. If a small force Fext > is 
applied, the filament velocity remains an approximately 
linear function of the applied force, but with an effective 



FIG. 6: The motor-induced sliding velocity Vs of an actin 
filament in the absence of external force is shown as a func- 
tion of e = io/Ss for N = 10 (dotted line), iV = 25 (dashed 
line), N = 100 (dashed-dotted line) and iV = 500 (solid 
line). We observe that Vs — >■ vo for stiff motors as A*' is in- 



creased. Parameter values: 
a/k = 1.33 pN. 



C = 0.002 pN (nm) 



0.06, 



friction greatly enhanced by motor binding/unbinding. 
This enhancement of friction is also termed in the lit- 
erature as protein friction [28]. At high velocity, only 
a few motors are attached to the filament and the fila- 
ment velocity approaches the value it would have in the 
absence of motors as the applied force is increased be- 
yond a characteristic value. When the external force is 
ramped down the filament velocity jumps to the lower 
branch corresponding to a lower value of the force, re- 
sulting in hysteresis. 



B. Fluctuation Dynamics 

We now examine the dynamics of noise-induced fluc- 
tuations about the steady state by letting Si — x ~ v, 
where v is the steady state velocity, given by the solu- 
tion of Eq. ( 23 ) discussed in the previous section. The 



dynamics of the fluctuation 5x is then described by the 
equation 



C(5i = ~kNA{v)Spb - kNpbsSA + ■q{t) 



(30) 



where both (5A — [dvA{v)]5x and Spb{t) depend on noise 
only implicitly through the velocity i, with 



dtSpb = 



1 

t{v) 



Spb -Pbs{v) 



d_ 
dv 



1 

t{v) 



Sx (31) 



The random force ri{t) in Eq. (30) describes noise on the 
filament, with (r;(i)) = a.nd{r]it)T]{t')) = 2BS{t - t'). 
Noise can arise in the system from a variety of sources, 
including the fluid through which the filament moves and 
the motor on/off dynamics. For simplicity we assume the 
spectrum is white, albeit with a non-thermal strength B. 
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By solving Eq. (31 ) with initial condition 5pi,{t = 0) = 



and substituting in Eq. ( 30 ) , we obtain a single equation 
for Si, 

[C + Ca{v)] Sx{t) + LOoCaiv) f dt' e-''^'-'"^ 5x{t') = r^{t) 

Jo 

(32) 

where we have introduced an effective frequency n{v) = 
T^^{v) + u!b and active frictions 



Ca(w) = kNpbs{v)dyA{v) 

d fl 



Civ) = kNpUv)A{v) 



dv 



(33) 
(34) 



In all the parameters defined above v has to be replaced 
by the steady state solution obtained in the previous sec- 



tion. The time scale 
cycle of a loaded motor 
Ca given by Eq. (pTl). It 



^ represent the duration of the 
Note that Ca{v = 0) = Ca, with 
is evident from Eq. (32) that 



motor dynamics yields a non-Markovian contribution to 
the friction. 



If we neglect the load dependence of the unbinding 
rate by letting = 0, hence r^^ — ujq, then Co(^^) = 
C„o = Nrki/vo and Ca(w) = 0. In this hmit {[5x{t) - 
(5x(0)]^) — 2Daot and is diffusive at all times, with an 
effective diffusion constant Dao 



B 



When v is finite we obtain 



[Sx{t) - Sx{0)f) = 2Dat + ADa 



1 2 



r 



(35) 



where Da = B/[C + Ca{v)]^ and flaiv) = n{v) + 
'^oCa(^)/[C + Ca{v)]- The characteristic time scale 
controls the crossover from ballistic behavior for t ^ fl"^ 
to diffusive behavior for t ^ ^a^- 1^ determined by the 
smaller of two time scales: ft'-^, defined after Eq. (32 1, 



that represents the duration of the cycle of a loaded mo- 
tor, and the active time (t^oC/IC + Ca])"^ that represents 
the correlation time for the effect of motor on/off dy- 
namics on the filament. At long times the mean-square 
displacement is always diffusive, with an effective diffu- 
sion constant 



Dr, 



[C + Ca{v)]na 



(36) 



This result only describes the behavior of the system in 
the stable region, where the effective friction remains 
positive. At the onset of negative friction instability 
C + Ca{v) —T' and the effective diffusivity diverges. In 
other words the instability is also associated with large 
fluctuations in he rod's displacements due to the cooper- 
ative motor dynamics. 

To leading order in u the frequency fla that controls 
the crossover to diffusive behavior is simply O ~ -I- 
ujt, + 0{v'^). For non-processive motors such as myosins 
loq ^ Wf, and ^l ^ uq. The effective diffusion constant is 
given by 



cff - 



2C CaO ( VqU \ . r, , ^z| 
(C + Cao)^ \UJo[l + e) 



(37) 

This expression indicates that the enhancement of the 
diffusion constant comes from the competition of the 



ballistic motor-driven motion of the filament at speed 
^ ^^oCao/(C + Cao) and the randomization of such motion 



by the motor on/off dynamics on time scales ■ 



The 



result is that the filament dynamics is diffusive at long 
times, but with an enhanced diffusion constant. 

Finally, we stress that the correlation function ([(5a;(i) — 
(5x(0)]^) describes the fluctuations about the steady state 
value vt. if we write x{t) — vt + 6x{t) the actual mean 
square displacement of the center of mass of the rod is 
given by {{x{t) - x{0))^) = vH^ + {[Sx{t)-Sx{0)]^) and is 
ballistic at long times in one dimension due to the mean 
motion of the rod. In addition, due to nonlinearity of 
the Langevin equation (17) the mean value (x) in the 



presence of noise will in general differ from the steady 
state solution vt obtained in the absence of noise due 
to renormalization by fluctuations {Fa{x,t)) — Fa{v,t). 
These fluctuations are neglected in mean field theory. 



IV. ACTIVE FILAMENT DYNAMICS IN TWO 
DIMENSIONS 

In two dimensions the coupled translational and ro- 
tational dynamics of of the filament is described by 
Eqs. (la) and (lb). It is convenient to write the instan- 



taneous velocity of the center of the filament in terms of 
components longitudinal and transverse to the long axis 
of the filament, r — -I- Vj^fi. Similarly the stretch is 
written as A = A||U + Aj^n, where (see Eq. ( |12b[ )) 

A|l = u- A = (V^l +t;„)T , (38a) 
A_L = fi • A = (VI -(- so^)r . (38b) 
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It is then clear that Ay has the same form as in one 
dimension 



All = 



(Vh - vq)/ujo 



(39) 



and the mean-field value of the attachment time t is given 
by 

(40) 

where we have carried out the average over sq. Inserting 
these expressions in Eqs. (15) and (16), the mean field 



active force and torque exerted by bound motors on the 
filament can then be written as 



r„ = -kNpbit) 

Ta = -kNpb{t)T 



(Vh - Vq)/u!q 



f-1 + e 



u + V±Th 



V±VmT 



(41a) 
(41b) 



Ca given in Eq. (27) for a rod in one dimension, with 
A — > A||. The transverse and rotational friction coeffi- 
cients are enhanced by motor activity. Their active com- 
ponents are given by 



C-La 



kNrTQ 



r + (1 - r)fo 1 



kNrToL^/12 
r + (1 - r)ro 



_ fcjVrroi;o(To-|-e|A^|) 



(43a) 
(43b) 

When the load 



Finally we have, ga 

dependence of the unbinding rate is neglected {u — 0), 
all friction coefficients are enhanced by motor activ- 
ity. When the force/ velocity and torque/angular veloc- 
ity curves are calculated to nonlinear order, we find that 
the only instability is the negative longitudinal friction 
instability obtained in one dimension. No instabilities 
are obtained in the angular dynamics. We expect this 
will change if we include the semiflexibility of the fila- 
ment 1291 [301. 



A. Steady State and its stability 



B. Fluctuations around the steady state 



The steady state of the motor-driven filament in two 
dimensions in the absence of noise is characterized by 
the center of mass velocity v = t;||U + uj^n and angular 
velocity i). In the absence of any external force or torque, 
"i? and v±_ are identically zero, whereas the longitudinal 
dynamics described by v\\ is identical to that obtained in 
one-dimension: the filament will slide along its long axis 
at a steady longitudinal velocity i'n = Fp/{C + Ca), with 
Fp and (a given by Eqs. (26) and (27), respectively. 

To gain some insight into the stability of the system 
under application of external forces or torques, we ex- 
pand Fa and to linear order in velocities v and as, 

Fa(v,^)^F,+ (f^)^.|, + (f^)^., + (^)^^,and 

^^.(v,^)- + + where 

Fp = Fa,o = -PpU, is the tangential propulsion force due 
to the motors. The subscript '0' indicates that the ex- 
pressions are evaluated at v = and -d = 0. This leads to 
steady state force/velocity and torque/velocity relations 
of the form 



iCe +C0a)^ = Text - 9aV± 



FpU 



(42a) 
(42b) 



where we have introduced an active "momentum" ga 
given by ga = — (f^) • The active contributions to 
the longitudinal, transverse and rotational friction co- 
efficients are defined as C||a — ~u ' (^jp) ' Cj-a — 

-n • (iff)^, and Cea = ~ (ff The longitudinal 
friction coefficient Ola is identical to the active friction 



We now examine the dynamics of noise-induced fluc- 
tuations about the steady state by letting bv ~ r — v 
and 56 — 9 — -d where v and d are the steady state ve- 
locity and angular frequency in the absence of external 
force and torque. As noted in the previous section when 
Foxt — and Text = 0, wn = w ^ 0, with v given by the 



solution of Eq. (23), and tjj^ = = 0. Projecting velocity 



fluctuations longitudinal and transverse to the filament, 
(5r = U(5Vj| -I- n5Vj_, the dynamics of fluctuations is de- 
scribed by the coupled equations, 

[C|l + C\\a[v)\ 5V\\ - -fciVAii {v)5pk{t) + r;|| , (44a) 
[a + aa{v)]5V^=r^^, (44b) 

[Ce + Cea(w)] 59 = -kNpbs{v)T{v)vra{,v)5VA^ + rie , 

(44c) 



with 



[Ce + Cea{v)] 5pb = ~Q.{v)5pb - Pbs{v) 



d_ 
dv 



1 



t{v) 



5Vn 



(45) 

where the effective frequency n{v) = t ^(v) + cub 
and the longitudinal active friction C\\a{v) are as in 
one dimension, C±a{v) = kNpbs{v)T{v) and Cea{v) — 
kNpbs{v)T{v)L'^ /12. In all the parameters, v = U|| has to 
be replaced by the steady state solution obtained in one 
dimension in the absence of external force or torque. 

The time-correlation function of orientational fluctu- 
ations, A9{t) = 59{t) — 59{Q), can be calculated from 
Eqs. (|44b| and (|44c|, with the result 



{M{t)M{t')) ^2Dea mm{t,t') 



(46) 



11 



The effective rotational diffusion constant is enhanced by 
the transverse diffusivity and is given by 

Dea[v) = ^ H ^ (47) 

[Ce + Cea{v)f [a+CUv)f 

with 4 (-d) = [Ce + Cea{v)] /kNpbs{v)T{v)vjn{v). Using 
Eq. ( 46 ) , one immediately obtains the angular time- 



correlation function as |31| . 



(48) 



The fluctuations in the probability of bound motors are 
driven by their coupling to the stochastic longitudinal 
dynamics of the filament. Assuming 5pf,{0) = 0, we ob- 
tain 



nr. 



-n„|t-t'| 



g-n„(t+t')l ^ 
(49) 



where ^ai^) = fl{v) + luq-^ 

Civ) = kNpUv)A^\{v)^{l), and v}v)' = 
iVfcA||(u)/ + C||a(^^)] is a longitudinal propulsion 
velocity. Notice that Vp{v = 0) = Vs/pbso, with Vs given 
in Eq. ( p8| ). Finally, we can compute the correlation 
function of the fluctuation 6r of the filament's position. 
In the laboratory frame the dynamics of Sr can be recast 
in the form of a simple equation, 



6r ^ ^VpSpb{t)u+ C + C{v) 



(50) 



Fluctuations in the probability of bound motors do not 
couple to orientational fluctuations to linear order. It is 
then straightforward to calculate the correlation function 
of displacement fluctuations, with the result 



(Mt) - Sr{0)f 



2D^s t + 



r 



41^2 



0a 



(51) 



where effective longitudinal and transverse diffusion con- 
stants have been defined as 



^||a=B||/[C||+C||a(l^)]', 



(52a) 
(52b) 



Finally, using r{t) = Sr{t)+ /J dt'vu{t'), the mean square 
displacement (MSD) can be written as, 



([r(t)-r(0)]2) = {Mt)-6vm 



V 



2 r 



-Deat 



Dea 



(53) 

The MSD is controlled by the interplay of two time scales, 
the rotational diffusion time, Dg^, that is decreased by 
activity as compared to its bare value, Dg^ , and the time 
scale which is turn controlled by the duration of the 
motor binding/unbinding cycle. If Dg^ ^ which is 
indeed the case for actomyosin systems [35] then on times 
t:$>n^^ the MSD is given by 



([r(t)-r(0)]2) 
with 

-Deft Diu + D±a 



1 - e 



Do 



(54) 



Dea + na V[C|| +C||a(")]f^a 



(55) 

In other words the rod performs a persistent random walk 
consisting of ballistic segments at speed v randomized 



by rotational diffusion. The behavior is diffusive both 
at short and long times, albeit with different diffusion 
constants, D^s and I?off + v'^ / {"^Dga), respectively. This 
is indeed the dynamics of a self-propelled rod. If the noise 
strengths B\\, B± and Bg are negligible, then Eq. (54) 
reduces to 



([r(t)-r(0)]^ 



V 



0a 



t 



1 



-Deat 



0a 



(56) 



and the MSD exhibits a crossover from ballistic behavior 
for t <C Dg^ to diffusive at long times. 

It is worthwhile to note that if one neglects load 
dependence of unbinding rate by taking v = 0, ef- 
fective diffusivity at long time is enhanced with, 
Dls = -Dya + + (^^°)V2^ea> due to the interplay 
between ballistic motion driven by the tethered motors 
and rotational diffusion, unlike the situation in one 
dimension. 



V. SUMMARY AND OUTLOOK 

We have investigated the dynamics of a single cy- 
toskeletal filament modeled as a rigid rod interacting with 
tethered motor proteins in a motility assay in two dimen- 
sions. Motor activity yields both an effective propulsion 
of the filament along its long axis and a renormalization 
of all friction coefficients. The longitudinal friction can 
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change sign leading to an instability in the filament's re- 
sponse to external force, as demonstrated by previous 
authors [TT]. The effective propulsion force and filament 
velocity in the steady state are calculated in terms of 
microscopic motor and filament parameters. 

We also considered the fluctuations of the filament dis- 
placement about its steady state value and demonstrated 
that the coupling to the binding/unbinding dynamics of 
the the motors yields non-Markovian fluctuations and en- 
hanced diffusion. Future work will include the stochas- 
ticity in the motor displacements and the semiflcxibility 
of filaments, which is expected to lead to buckling insta- 
bilities [32] and anomalous fluctuations |33| . 
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Appendix A: Solution of Mean-Field Equation 



Here we discuss the solution of the mean-field equation 



( 13 1 for the attachment time t. For simplicity, we con- 
sider the one-dimensional case in detail. The discussion 
is then easily generalized to two dimensions. The mean- 
field equation for the residence time r is rewritten here 
for clarity: 



tmf 



w-1(A(tmf)) 



(Al) 



The solution clearly depends on the form chosen to de- 
scribe the dependence of the motor unbinding rate on the 
stretch A, in turn given by /S.{tmf) = (i — vq)/{t^Ip + 
ewo). The mean-field equation must be inverted to de- 
termine Tmf as a function of the filament velocity x — v. 
For compactness we drop the label 'MF'. It is clear that 
T has a maximum at f = wo, where t = Wq"^- This sim- 
ply corresponds to the fact that the time a motor protein 
spends attached to the actin filament is largest when the 
motors' tails are unstretched (A = 0) and the motors 
advance at the unloaded motor velocity, wq. 

It is convenient to use the dimensionless variable and 
parameters introduced in the text and write the stretch 
A as 



A = 



{u - 1)4 



1 



(A2) 



where u = v/vq, £D„ = uJu/coq and Iq = Vo/ujq. A form 
commonly used in the literature is the exponential form 



FIG. 7: Mean field attachment time tmf as a function of v for 
parameter values appropriate for acto-myosin systems: vo = 
1000 nm s~\ k = 10 pN nm~\ fs ^ 4 pN, a"^ = 7.5 nm, 
ojo — 0.5 (ms)~^, r = 0.06, corresponding to e = 5. The 
dashed line is the numerical solution of Eq. ( |A1[ ) obtained 
using the exponential dependence of the unbinding rate on 
the stretch. The solid line is obtained using the paraboUc 
ansatz given in Eq. ( A3 1 . 



a;„(A) = wqc"''^!, with a^^ a characteristic length scale. 
The dimensionless combination aA can then be written 
in terms of the parameter v — ai — a4/(l 4- e) and 
setting V = corresponds to neglecting the load depen- 
dence of the unbinding rate. The numerical solution of 
Eq. (Al) for the mean attachment time as a function of 
V is shown as a dashed line in Fig. [7] for parameter values 
appropriate for acto-myosin systems. As expected it has 
a sharp maximum at v = vq. At large v the attachment 
time decays logarithmically with velocity. As a result, the 
stretch is found to saturate at large velocity, as shown by 
the dashed curve in Fig. [8j This behavior is unphysical 
as it does not incorporate the fact that when the stretch 
exceeds a characteristic value of the order fd/k, the mo- 
tor head simply detaches, as shown in Fig. [2] Instead 
of incorporating this cutoff by hand, we have chosen to 
use a simple quadratic form for the dependence of the 
unbinding rate on the stretch, given by 



(A3) 



With this form the mean field equation (|Al| can be 
solved analytically, although the explicit solution is not 
terribly informative and will not be given here. The re- 
sulting attachment time is shown as a solid line in Fig. [7] 
The quadratic form reproduces the sharp maximum of t 
at V — Vq and yields r ~ at large v. The stretch 

then decays with velocity, as shown in Fig. [Sj 
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FIG. 8: Stretch A as a function of velocity v obtained us- 
ing the mean-field value of the attachment time displayed in 
Fig. [T] The parameter values are the same as in Fig. [7] The 
dashed line is obtained using the exponential dependence of 
the unbinding rate on the stretch. The solid line is obtained 
using the parabolic ansatz given in Eq. (A3 1. 
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